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A.  THEORETICAL  DEVELOPMENT  OF  NOVEL  OPTICAL  NIM 


I.  INTRODUCTION 

Based  on  modern  advanced  methods  of  Quantum  Optics,  we  discuss  here  the  possibility 
of  achieving  the  negative  index  of  refraction  in  a  semiconductor  with  donor-like 
impurities.  This  approach  is  genuine  microscopic,  and  uses,  instead  of  Maxwell  equations 
(homogeneous  NIM  designs),  the  Quantum  Mechanical  methods.  This  design  requires 
several  specially  prepared  laser  beams,  one  probe  beam  and  several  support  beams,  with 
specific  intensities,  frequencies  and  polarizations.  These  laser  beams  pump  energy  levels 
and  create  the  desirable  density  matrix  elements  of  an  interest,  which,  in  their  turn, 
produce  coherent  penneability  and  permittivity  responses  at  frequency  of  the  probe  laser 
beam.  The  quantum  states  of  hydrogen-like  donor  atom  embedded  into  solid  state  matrix 
and  states  of  an  electron  in  conduction  band  constitute  a  discrete-level  atomic  medium 
within  the  optical  range.  The  coherent  coupling  of  an  electric  dipole  transition  with  a 
magnetic  dipole  transition  leads  to  resonance-like  penneability  and  pennittivity 
responses,  which,  within  some  frequency  band,  ensures  the  negative  refractive  index.  The 
implementation  of  this  scheme  is  carried  out  in  Sn-Zn-doped  indium  oxide,  In2-xSnx03 
:Zn  (Zn-doped  ITO),  and  the  calculations  show  feasibility  of  the  effect  within  a  broad 
bandwidth  with  an  extremely  high  figure  of  merit  ( FOM  =  |  Re(«)/hn(«)  |  >10)  in  the 
THz  regime. 


The  main  ideas  behind  the  Quantum  Optics  approach  which  have  developed  recently  will 
be  described  next.  Several  years  ago,  Oktel  and  Mustecaplioglu  [1]  suggested  an  elegant 
method  to  achieve  the  negative  refractive  index  effect  based  on  Quantum  Optics.  Despite 
the  very  different  physics,  their  method  is  similar  to  electromagnetically-induced 
transparency  formalism  (EIT)  [2],  Their  model  suggests  the  following  arguments.  The 
main  difficulty  of  reaching  the  negative  refractive  index  in  the  optical  regime,  in  a 
homogeneous  medium,  is  the  frequency-dependent  permeability  function,  //(<»),  which 
in  non-ferromagnetics,  is  close  to  1  with  no  magnetic  response.  This  behavior  was 
explained  by  Landau  [3],  who  showed  that,  in  the  optical  regime,  the  projection  of  the 
magnetic  moment,  fi  ,  of  an  individual  spin,  s,  onto  the  high-frequency  magnetic  field, 

H(oo),  is  almost  zero  due  to  the  low  frequency  of  the  spin  precession,  a>s  =  2sjubH  I  h 
( juB  is  the  Bohr  magneton),  compared  to  the  frequency  of  the  e.m.  wave,  co .  Note,  that  in 
non-homogeneous  situation,  this  difficulty  is  avoided  (at  the  expense  of  high  scattering 
losses)  due  to  the  effective  nature  of  /j{co)  in  such  a  way  that  the  displacement  currents 
or  Mie  resonances  cause  the  considerable  magnetic  response  which  is  effectively 
equivalent  to  //  ^  1 . 

We  report  here  solid-state-based  negative  refractive  index,  homogeneous,  and  low-loss 
material.  Obviously,  as  was  mentioned  in  Ref.  [1],  since  N  in  the  solid  state  can  be  ~  3-4 
orders  of  magnitude  larger  than  N  in  the  gas  phase,  it  would  be  extremely  desirable  to 
implement  the  solid  state  design.  The  solid  state  implementation,  in  principle,  can  allow 
for  the  desirable  negative  refractive  index  effect  at  THz  (and  higher)  frequencies  (the  gas 
implementation  [4-5,  6-9]  so  far  provides  the  effect  within  MHz  range,  and  not  at  THz 
range).  Also,  since  the  density  of  transitions  is  much  higher  in  the  solid  state  than  in  the 
gas  phase,  one  might  anticipate  better  parameters  of  the  negative  refractive  index  band  in 
tenns  of  the  magnitude  of  the  refractive  index  and  the  bandwidth. 

This  solid-state-based  negative  refractive  index  material  is  based  on  a  n-type 
semiconductor  doped  simultaneously  with  a  transition-metal  (TM)  with  an  atomic 
volume  density,  N,  and  p-type  acceptor,  with  approximately  the  same  density,  N.  The 
desirable  negative  refractive  index  effect  is  electromagnetically  induced  by  the  adjusted 
support  e.m.  fields  (laser  beams)  to  ensure  the  coherence  of  the  dielectric  and  the 
magnetic  responses.  In  accordance  with  band  theory  [10-12],  the  donor  atom,  or  “defect”, 
produces  the  discrete  energy  levels  within  the  gap  between  the  valence  band  and  the 
conduction  band  (CB)  (Fig.  1).  The  appropriate  states  of  a  donor  can  be  treated  as  states 
of  the  hydrogen-like  atom  with  re-nonnalized  effective  mass,  m  ,  of  an  electron  and  re¬ 
normalized  effective  charge,  e  ,  of  an  electron.  Three  first  lowest-energy  donor  states  are 
supplemented  by  the  additional  4-th  electronic  state  in  the  conduction  band,  to  couple  the 
donor  states  with  the  states  of  a  crystal.  The  required  non-zero  magnetic  moment  between 
two  hydrogen-like  atom  states  with  the  same  parities,  but  with  mutually  orthogonal  radial 
wave  functions  is  achieved  by  the  specially  oriented  additional  support  beam  which  acts 
as  a  perturbation.  Consequently,  we  explore  the  opportunity,  by  adjusting  the  parameters 
of  the  non-linear  system  (i.e.  the  density  of  the  dopants,  N,  the  intensities  of  e.m.  fields  of 


the  lasers,  etc.)  to  achieve  the  desirable  effect  in  the  THz  domain.  All  of  the  parameters, 
as  we  have  mentioned  above,  are  implemented  for  the  specific  example  of  tin-Zn  doped 
indium  oxide,  In2-xSnx03  :Zn  (Zn-doped  ITO),  to  demonstrate  the  feasibility  of  the 
suggested  design.  In  a  detailed  theoretical  analysis,  we  have  derived  the  quantum 
Louiville  equations  for  the  density  matrix  elements,  ptj ,  of  the  atomic  states  (indexes  i 

and  j  denote  atomic  levels).  Based  on  methods  of  Quantum  Optics,  these  equations  were 
solved  and  the  solutions,  for  non-diagonal  matrix  elements  with  appropriate  coherences 
are  analyzed.  Next,  we  have  calculated  the  relevant  dipole  moments  for  the  allowed 
dipole  transitions  between  states  with  different  parities.  Afterwards,  based  on 
perturbation  theory,  we  have  calculated  the  induced  non-diagonal  magnetic  moment 
element  responsible  for  the  transitions  between  states  with  similar  parities.  All  of  these 
density  matrix  elements,  and  the  dipole  and  magnetic  moments  are  utilized  to  calculate 
the  dielectric  response,  £{co) ,  the  magnetic  response,  ju{co) ,  and  the  refractive  index, 
n(eo) . 


II.  Atomic  level  system 

In  accordance  with  the  band  theory  [10-12],  in  a  n-type  doped  semiconductor,  if  the 
density  of  dopant  atoms,  N,  is  not  too  large,  the  states  of  the  atoms  are  mutually 
independent  and  they  correspond  to  well-isolated  discrete  energy  levels  within  the  band 
gap.  The  appropriate  first  3  energy  levels  are  shown  in  Fig.  1.  The  dopant  atom,  or  zero- 
order  defect,  donates  electrons  of  its  outer  electron  shell  into  the  conduction  band  (CB), 
and  consequently,  is  effectively  charged  with  some  integer  charge  Z  (Z  =1,  2  or  3), 
depending  on  the  type  of  the  crystal  and  the  type  of  a  defect  (dopant). 


Fig.  1.  System  with  discrete  energy  levels  as  described  in  text. 


The  long-range  behavior  of  the  defect  potential  V(r)  is  described  by  the  Coulomb-like 
potential  of  the  electron  with  re-normalized  effective  mass,  m0  — >  in  ,  and  re-normalized 
electron  charge,  e  — >  el  ,  (e0  is  a  static  dielectric  constant)  [10-11]  .  The  appropriate 
Schrodinger  equation  provides  the  hydrogen-like  wave  functions,  the  Bohr  radius, 


ah  =  fr2s0/  Zrne1 ,  and  the  energy  levels,  E  =  -Z2m*e4  /  2h2s2n2 ,  with  nq  =  1,2,3,....  as 
main  quantum  numbers.  We  are  interested  in  the  first  3  lowest  energy  states  ( nq  =1,2,3) 
which,  in  the  text,  are  called  0-state  (ground  state,  nq  =1),  1 -state  (second  excited  state 
nq  =  3 ),  and  2-state  (first  excited  state  nq  =  2  ).  All  energies,  E. ,  of  these  levels  are 
below  the  bottom  of  the  conduction  band  as  shown  in  Fig.  1.  The  orbital  quantum 
numbers,  /,. ,  of  these  states,  which  define  the  parity,  P.  =  (-l)/f  ,  are  supposed  to  be: 
l0  =  0,  /j  =  l2  =  1 ,  to  ensure  the  same  parity  of  1  and  2  states,  which  is  different  from  the 
parity  of  the  0-state.  The  4-th  energy  level  (3-state)  of  our  system  is  an  electron  energy 
level  separated  by  the  same  distance  from  level  1  as  level  2  (i.e.  E3  -  E]  =  El  -  E2 ),  and  it 
has  the  same  parity  as  the  0-state  with  l3  =  0  .  The  magnetic  quantum  numbers  of  all 

four  states  are  the  same  (m  =  0).  The  3 -state  is  inside  the  conduction  band,  and  this  state 
connects  the  isolated  hydrogen-like  atom  states  with  the  CB  states  of  a  solid-state  matrix. 
The  following  electric  dipole  operators  are  non-zero: 

<  0  |  d  1 1  0,  <  0  |  d  |  2  0,  <  2  |  d  |  3  0  ,  in  accordance  with  the  parities.  The 

magnetic  dipole  operator  <1|//|2>^0  between  1 -state  and  2-state  can  be  non-zero, 
since  parity  rules,  but  it  is  actually  zero  in  the  undisturbed  situation  due  to  the 
orthogonality  of  the  radial  wave  functions.  Next,  we  introduce  the  four  external 
electromagnetic  fields  which  are  to  couple  with  these  states.  We  assume  that  the  probe 
e.m.  field  (beam  b2 )  provides  a  strong  electromagnetic  coupling  for  the  electric  dipolar 
transition  |  0  >  — >  |  2  >  (the  standard  symbol  |n>  stands  for  the  state  n).  This  field  is 
described  by  the  amplitude,  E02 ,  of  the  electric  field  with  frequency,  co02 ,  detuning,  S02 , 
and  the  Rabi  frequency,  Q02  =  d02E02  / ti .  The  support  e.m.  field  (beam  bx)  is  responsible 
for  the  weak  electromagnetic  coupling,  Q01  =  dmE0X  /  h  ,  due  to  electric  dipolar  transition 
1 0  > — >|  1  >  with  frequency,  co01,  detuning,  A,  and  the  Rabi  frequency,  Q01  =dmE0l/h 
(l«oi  l«l  Q02  | ).  The  driving  e.m.  field  (beam  b3 )  provides  the  weak  electromagnetic 
coupling  due  to  electric  dipolar  transition  |  2  >  — >  |  3  >  between  the  atomic  state  and  the 
solid-state  matrix,  with  frequency,  co23  ,  detuning,  S2i  ,  and  the  Rabi  frequency, 

Q23  =  d23E23  /  h  .  All  these  e.m.  fields  are  assumed  to  posses  the  same  polarizations  with 
collinear  Poynting  vectors,  501  \\S02  ||523,  and  mutually  parallel  electric  and  magnetic 
fields  (E0l  ||  E02 1|  E23 ,  Hm  ||  H02 1|  Pl2 3). 


III.  Density  Matrix 

Such  a  system  is  described  by  a  Hamiltonian  in  the  fonn  [5]: 

H  =  H0  +  Hi  0) 

where  H0  is  a  field-free  Hamiltonian  and  Hi  is  the  interaction  with  the  e.m.  fields: 


(2) 


(3) 


<k\H,  \l>  =  --nnkIe-K,0t-(i>,), 

'  2 

Here  the  Rabi  frequencies  associated  with  optical  transitions  are  defined  as 
Clkl  =  dkl-E/fi  where  dkl  =  e<k\r\l>  is  the  electric  dipole  operator,  e  is  the  charge  of 

the  electron,  and  E  is  the  amplitude  of  the  relevant  electric  field.  In  the  rotating- wave 
approximation,  within  the  semiclassical  theory  of  optical  interactions  [1,  5],  the  density 
matrix,  p  ,  of  the  system  evolves  as  required  by  the  Louiville  equations: 


^-=-Uh,p]-Ut,p) 

dt  n  2 

(4) 

(5) 

k,l 


with  diagonal  elements  of  the  relaxation  matrix  <  i  \  F  \  i  >  =  yu  and  off-diagonal  elements 
defined  as  <  i  \  T  \  k  >  =  (ya  +  ykk)  / 2  ;  the  standard  simbols  [A,  B]  and  { A,B}  stand  for  the 
commutator,  and  symmetrized  product  of  operators  A  and  B,  respectively.  The  detunings 
are  defined  by  tiSkl  =  E,  -Ek  - hcokl  with  cok,  as  the  angular  frequency  of  the  appropriate 
e.m.  field,  and  the  detuning  of  the  support  field  is  called  8m  =  A ,  so  there  are  a  total  of 
three  detunings:  S02 ,  S23  and  A . 

Next,  we  need  to  obtain  the  solutions  for  p20,  pl0,  pl2,  p23,  and  pn  of  Eqs.  (4), 

which  affect  the  electric  and  the  magnetic  responses,  and  we  will  follow  the  EIT 
approximation  [2]  which  is  based  on  the  usage  of  small  parameters: 
|Q0I/Q02|«1;  |Q23/Q01|«1  .  The  equations  for  p20,  pw,  pX2  are  split  off  from 
the  rest  of  Eqs.  (4)  due  to  the  strong  interactions  between  0,  1,  and  2  levels.  The  system 
of  three  equations  was  solved,  and  the  solutions  for  pXQ ,  p20  and  pl2  are  given  by: 


A  o 


A  2 


_  f202  exp[  i(ojl}2  S02  )t] 

(6) 

(A-S32-iy22) 

zQoi[ri2  +  Z'(A  -  £02)]  exp[-/(®01  -  A)r] 

(7) 

[(Tii  +z'A)(y12  +z(A-^02))+  Q02  2] 

Q01Qq2  exp[-z'(A<»  +  S02  -  A)t] 

(8) 

[( Yu  +Z'A)(/l2  +z'(A-cS’o2))+  |  f^02  I"] 


rii  =  (An  +  r00)/2;  r'n  =  (r22  +r00)/2i  r\2  =  (Ai i  +  a22 ) / 2 . ^02  =  (e2-e0)/^ 


(9) 


and  A a>  =  (El-E2)/h  is  the  energy  difference  between  level  1  and  level  2.  These 
solutions  lead  to  the  following  non-diagonal  element  of  the  density  matrix,  p, , : 


(10) 


O*01Q02Q23  exp[-z'(A<»  +  A- S02  -S2i)t~\ 

Pa  (kco-iyn)[(yu  +  iA)(yn  +/(A-A02»+ 1 Q02 12] 

Since  we  need  to  obtain  an  electric  dipole  and  a  magnetic  dipole  oscillating  at  the  same 
frequency  in  order  to  modify  the  effective  pennittivity  and  permeability  responses 
simultaneously,  the  coherence  conditions  require  synchronization  of  frequency-dependent 
factors  in  Eqs.  (6,  8,  10): 

co02  -  ^02  =  A®  +  S02  -  A  =  A co  +  A  -  S01  -  S22  (11) 

Equation  (11)  insists  that  the  tunings  are  actually  not  independent: 

A  =  Aco  -  co02  +  2A02  (12) 

S2>=2(A-St2)  (13) 

Equations  (12)-(13)  express  tunings  A  and  S2i  as  functions  of  the  tuning  of  the  probe 
beam,  S02 .  Eqs.  (12-13)  correspond  to  the  resonance  situation  when  the  e.m.  field  of  the 
probe  beam  has  frequency  O)prohe  =  3Z2/n  e4  /8/F.y]  -S02  and  the  e.m.  field  frequency  for 
0-1  transition  is  equal  to  the  frequency  for  2-3  transition,  com  =  oo23  =2 coprobe,  and  it  is 
twice  larger  than  the  co  be .  Hence,  all  four  of  the  effective  energy  levels  are  equidistant, 

and  the  support  beams  and  the  driving  beams  are  degenerated  since  they  have  the  same 
frequency.  As  a  result,  we  will  exclude  the  driving  beam  and  leave  only  support  beam 
with  frequency  2 co  ^  with  substitution  Q23  =Q01  into  Eq.  10.  Also,  the  probe  beam 

causes  two  additional  transitions,  one  is  electric  dipolar  between  levels  1  and  3,  and 
another-  magnetic  dipolar  transition  between  levels  2  and  1 .  The  single  tuning,  S02 , 
defines  the  frequency-dependence  of  both  electric  and  magnetic  responses. 


IV.  Permittivity  and  permeability  functions 

The  classical  electric  dipole,  P ,  is  calculated  from  the  general  expression  [13]  as: 

P  =  Tr{pd)  (14) 

where  d  is  the  dipole  operator.  The  amplitude  of  the  electric  dipole  oscillating  at 
frequency  co  =  co02  of  the  probe  e.m.  field  is  defined  by  the  0-2  and  1-3  transitions: 

P  —  d02p20  +  dup2l  (15) 

The  effective  electric  polarizability  second  rank  tensor,  ae,  is  defined  by  its  Fourier 
transformation: 


(16) 


?{co)  =  ae(co)E(co ) 

where  E  is  the  electric  field  of  the  probe  beam.  Since  we  are  interested  in  a  single 
frequency  (co  =  co02 ),  we  will  drop  the  co  dependence,  and  let  ae{co)  =  ae .  Since  we  will 

consider  d  parallel  to  the  E  direction,  ae  is  a  scalar  calculated  from  Eqs.  (16),  (6)  and 
(10)  as: 

1  J  I  Q  |2  7  2  (17) 

a  _  _ a\3a02  I  ^01  I _ I _ ^02 _  V  ' 

h(A(0 - iyl3)[(y\ I  +  i^Kr'n  +  i(A- S02))+  |  Q02 12]  h(co02-Aco-iy22) 

We  should  mention  that  the  second  term  in  the  right  part  of  Eq.  (17)  does  not  depend  on 
the  tuning,  d02 ,  and  this  tenn  is  much  smaller  than  the  first  tenn.  We  also  have  included, 

into  the  expression  for  the  effective  pennittivity  s  =  1  +  \naN  ,  D  =  sE  (electric 
displacement),  the  Drude  background  permittivity,  which  is  due  to  the  free  carriers 
(electrons)  in  the  CB: 


s  = 


47T  Ndl3d02 


Cl 


01 


h(Aa)-iy  n)[(yn+iA)(yn  +i(A-S01))+  \  Q02 1‘] 
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co  +iy  co 


-) 


(18) 
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Here  the  Drude  plasmon  frequency,  co  =  - 2 —  ,  depends  on  the  concentration,  nc ,  of 


s  m. 


the  free  carriers  in  the  CB,  yp  ~0.lcop,  and  we  use,  in  our  calculations,  the  following 
concentration:  nc  ~  1 0 1 2  (1  /cm3)  . 


The  effective  permeability,  ju  ,  is  calculated  from  the  classical  magnetization,  M  : 


M  =  Tr{pm }  (19) 

where  m  is  the  magnetic  moment  operator.  Again,  we  are  interested  in  m  =  Mn  parallel 
to  the  magnetic  field  H  of  the  probe  beam,  and  pn  is  the  relevant  p  matrix  element  in 
Eq.(19).  The  effective  p ,  obtained  analogous  to  the  procedure  described  above,  is  given 
by: 


+  nM„Nimd’aE,t  (20) 

h2[(yn+iA)(y'l2+i(A-S02))+\Cl02 12] 

Here,  we  assume  a  standard  zero  background  magnetic  response  situation  since  no 
magnetic  response  to  the  field  H  within  the  optical  domain  is  possible  for  non¬ 
ferromagnetic  situation  [3]  . 


V.  Dipole  and  magnetic  moments 

To  calculate  the  dipole  moments,  dm,  d02,  dl3,  the  wave  functions  of  the  atomic  states  are 
needed.  These  wave  functions,  are  the  products  of  the  radial  wave  function  Rnl  (r)  and  the 
angular  wave  function  Ylm(0,cp )  [13]: 


^n,^r,d,(p)  =  Rnl{r)Ylm{e,(p) 


(21) 


The  relevant  angular  wave  functions  are:  Ym  =  1  /  Vdvr  (0-state,  3-state), 
Yw(6,<p)  =  i cos#V3/ 4 n  (1 -state,  2-state),  and  the  radial  wave  functions  are  given  by: 
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(22) 

(23) 


with  a  =  ab  and  p  =  2Z/7  na  .  The  calculation  provides  the  explicit  wave  functions: 


1  27  7r 

'J'ioo (r,e,<p)  =  -j=( — )3/2  exP( - )  (0-State) 

V8  n  a  a 


^210  ir,e,(p)  = 
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2V8 n  a  a  2 a 


j  9  ?  7  9  7r  9  7r  7r 

^3  w(r,e,cp)  =  -x  —  (— )3/2(— )(4-—  )exp(-—  )cos£  (1-state) 
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(24) 

(25) 

(26) 


The  electric  dipole  operator,  d=er  ,  has  a  dz=ercos0  component  in  z-direction 
(spherical  system  of  coordinates)  which  coincides  with  the  direction  of  the  electric  field, 

E ,  of  the  probe  beam.  Consequently,  the  dipole  moments  are  given  by  the  following 
integrals: 


2  n  K  oo 


d  02 


=  |  J  J vF*100(r, 6>, ^s>)er3  cosOW  2l0(r,  6,  (p)  sin  0d(pd0  dr 


0  0  0 

2n  n  oo 

d0l  =d*n  =  |  ||vl/*1o0(r,6,,^)er3  cos  6  VF310(/-,  6,cp)  sin  Od <pd 6 dr 
0  0  0 


(27) 

(28) 


The  dipole  moments  d0l  and  d02  calculated  from  Eqs.  (21-28)  are  below: 

,  ,  8 1  \j2iea, 

dm=d  13  = - 
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The  magnetic  dipole  transitions,  due  to  the  parities,  are  allowed  between  states  1  and  2  if 
the  magnetic  dipole  moment  between  these  states  is  non-zero  M21  0  (we  assume  the 

same  spin  wave  functions  Xm  f°r  all  states): 

2rrr  (31) 

M2i  =2Mb)  2w\r’d)*indr2d(Pdddr(<Xii2SX-U2) 

0  0  0 

with  the  spin  factor  Xin^Xm  =  1  /  2  . 

Since  states  1  and  2  belong  to  the  different  shells,  and  the  radial  parts  of  wave  functions, 
vF210(r,6l)  and  T/310(r,^)  ,  are  orthogonal  to  each  other,  the  r-integral  in  the  right  part  of 

Eq.(45),  and,  hence,  the  magnetic  moment,  M21 ,  both  are  equal  to  zero.  To  cope  with  this 
problem,  we  will  consider  the  induced  magnetic  moment,  M21  ,  due  to  the 
inhomogeneous  electric  field,  E ,  which  causes  the  dipole  moment-type  perturbation  V  : 

V  =  erE  (32) 


Here  rE  is  the  scalar  product  of  vectors.  This  field  is  oriented  perpendicular  to  the  other 
e.m.  fields,  in  such  a  way  that  its  electric  and  magnetic  fields,  E  ,H  ,  and  the  Poynting 
vector,  S  ,  are  mutually  perpendicular  to  E0l,H01  and  501 ,  appropriately.  We  assume 
that  the  relevant  Fourier  component  of  the  electric  field  has  an  amplitude,  E0,  and  the 
fields  E  =  E0  cxpiikr  +icot) ,  H  =  H0  cxp(ikr  +  icot )  are  at  a  low  frequency  compared  with 
the  probe  frequency  (i.e.  co«co0l).  In  the  following  text,  it  is  convenient  to  rotate  the 

initial  system  of  coordinates  (rotation  angle  0  =  n  /  2  )  in  such  way  that  the  new  z-axis  is 
perpendicular  to  the  old  z  quantization  axis,  and  the  appropriate  e.m.  field  of  the  probe 

and  the  perturbation  e.m.  field  are  perpendicular  to  each  other:  E  _L  Ep  ,H  _L  Hp  .  This 
transformation  allows  for  the  non-zero  magnetic  numbers  m  =  —1,  —0, +1  in  the  new 
system  of  coordinates,  and  the  perturbation  is  expressed  as: 

V  =  eE(r,  0)r  sin  9(ei,p  +  e /  2  (33) 

Next,  we  will  expand  the  long-wavelength  perturbation  into  the  series  ( kab  «1)  with 
accuracy  to  two  first  terms: 

E  «  E0(\  +  ikr  ~^(kr)2)  ^ 

We  assume  the  complex  wave  vector,  k  =  k  +ik  ,  and  the  relevant  part  of  the 
perturbation  which  contributes  into  the  matrix  element  M1X  is  shown  below: 


V  = -eE0(k  -ik  )r2  sin 9 cos 9(e,ip  +e  ,<3)/2 


(35) 


The  perturbation  above  leads  to  the  following  matrix  element,  M12 ,  between  state  1  and 
state  2  in  accordance  with  the  perturbation  theory  [13] : 

M  V '  V,(V*  +  Vil)M\  <36> 

21  Lu  E.-E- 

i,k 

with  three  diagonal  eigenvalues  M^  =  -1,0,+1  (magnetic  quantum  numbers)  and 
indexes,  i  and  k,  which  belong  to  1 -state  and  to  2-state,  respectively.  The  calculation  of 
Vik  from  Eq.  36  with  the  angular  wave  function  factors 

Yl0  =  i  cos  9^3/ 4k,  Yk±[  =  +/ V3  /  8tt  sin  9e±up  yields  the  magnetic  moment  of  an  interest: 

0.62  eahE0fuB{k' -ik)  (37) 
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VI.  Negative  refractive  index  band. 

The  complex  refractive  tuning-dependent  index,  n  =  n(S02 )  =  ±^£(S02)jU(S02  )  ,  with  the 
correct  sign  defined  by  the  conditions  [14-15]  : 

(|  s  |  -Re(^))(|  // 1  -Re(//))  >Im(£,)Im(//)  (38) 

is  expressed  as  a  function  of  the  probe  frequency  copmbe  =  <x>02  -  S02  (or  tuning  S02 )  with 
the  following  S02  -dependence  of  both  £  and  ju  : 
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The  resonance-like  function  F(802)  which  enters  Eqs.  (39-40)  is  given  by: 


F{S02)  = 


1 


[(y  a  +  i{Aco  co0 2  +  2S02)(y  2I  +  i{Aco  cd02  +  ^2))+  |  O02  |  ] 


(41) 


One  can  see  from  Eq.  (41)  that  F(S0 2)  has  two  resonances  at  the  frequencies  S02  =  8+es : 


3  1  1 

3 res  ~  —  _  ^02)  —  -y  _  ^02)  I  ^02 


(42) 
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The  behavior  of  the  function,  Re(  F(x )  )  with  dimensionless  parameter  x  =  S02  /  S„2  - 1  is 
shown  in  Fig.  2. 

F 


Fig.  2.  The  resonance  function  Re(F(x))  (Eq.  (41))  in  dimensionless  coordinate  of  probe  tuning 


One  can  see  that  Re  ( F(x)  )  has  two  resonances.  If  |  O02 1«|  A co-a>02 1,  the  right-wing 
solution  8+res  =ro02-Aty  =  llZ2m*e4/36^3£'o  ,  and  the  other  two  tunings  are: 
A  =  S^es ,  S32=  0.  This  resonance  of  interest  exactly  corresponds  to  the  frequency 
(orprobe  =  A®  =  5Z2m*e4  /  llfv’sl  of  the  probe  e.m.  field  equal  to  the  distance  between 

levels  1  and  2.  Expressions  (38)-(42)  were  utilized  to  calculate  both  the  real,  Re(«),  and 
imaginary  ,  Im (n),  parts  of  the  refractive  index,  n,  along  with  the  figure  of  merit, 
FOM=|Re(n)/Im(n)|,  as  a  function  of  the  probe  frequency  copmbe . 

The  specific  material  of  choice  is  the  transparent  conducting  tin-Zn-doped  indium  oxide, 
In2-xSnx03  :Zn  (ITO  doped  Zn)  [16-19],  with  the  following  relevant  intrinsic  parameters: 
Z=l,  <?0w4.4  [18]  ,  m*  &0.3me  [19]  ( me  is  the  mass  of  an  electron),  and  a  modest 

concentration  of  the  carriers  (electrons)  in  the  CB  nc  ~  10ls  (1/c/n3) ,  the  level  of  doping 
for  both  In  and  Zn  to  be  N  =  10 21  cm  3 .  The  level  of  losses  y  n  ~  yv  ~  108  Hz  was  taken 
as  lowest  in  gas  phase  [5],  and  then  it  was  increased  to  a  much  higher  level 
Zn  ~  Zii  ~  It)10  Hz,  to  check  the  stability  negative  refractive  index  effect.  The  free 
parameter  of  the  calculations  was  the  probe  beam  laser  intensity,  Ip ,  which  was  changing 
in  limits  of  experimentally  established  probe  laser  intensities 
10 6  W  / cm2  <  I  <  10s  W / cm 2  [20-21],  First,  we  calculated  the  frequency  dependence  of 
the  complex  permittivity  s{coprohe )  and  the  permeability  ^{coprobe)  which  are  shown  in 
Figs.  3-4. 


Fig.  3.  The  real  (blue),  Re[f(ffl  o^g)] ,  and  the  imaginary(red),  Im  [ff(®  o^g)] ,  parts  of  the  permittivity 

'  '  10 

as  a  function  of  the  probe  frequency  COprobe  .  This  is  high  loss  situation,  Y  Y\  2=10  Hz  ,  with  probe 
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laser  intensity  Ip  =  10  ( IF  /  cm  ) ,  and  doping  concentration  .V  =10  cm 
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Fig.  4.  The  real(blue),  Re[/r(©  ie)] ,  and  the  imaginary) red),  Im[//(®  feg)] ,  parts  of  the 

permeability  as  a  function  of  the  probe  frequency  COprobe  •  This  is  high  loss  situation, 

i  i  10  82 

y  n=  Yy  2  =10  /Tz  ,  with  probe  laser  intensity  Ip  =10  (W  /  cm  ) ,  and  doping  concentration 

21  -3 

N  =  10'  cm  . 

Fig.  5  demonstrates  the  frequency-dependent  complex  refractive  index  calculated  for  the 
probe  laser  intensities:  I  =10 8(W/cm2) ,  which  is  within  the  modern  limits  of  laser 
intensities.  The  negative  refractive  index  band  at  ~  12-13  THz  is  clearly  seen  with  a  well 


FOM 


pronounced  real  part  of  the  index  of  refraction  (Re(«) — 1)),  and  a  high  FOM  >  10  (Fig. 

6). 


Fig.  5.  The  real  (blue),  Re(«(<y)) ,  and  the  imaginary  (red),  Im(/?(®)) ,  parts  of  the  refractive  index  n  as 
function  of  the  probe  frequency  (Sn-Zn-doped  indium  oxide),  all  parameters  are  the  same  as  in  Figs.  3-4. 


Frequency  Hz 

Fig.  6.  Calculated  figure  of  merit  (FOM)  as  function  of  the  frequency  in  Sn-Zn-doped  indium  oxide,  all 

parameters  are  the  same  as  in  Figs.  3-4. 


B.  EXPERIMENTAL  DEVELOPMENT 


I.  Introduction 

The  experimental  approach  to  homogeneous  NIMs  included  two  different  material 
systems:  (Ini_xSnx)203  and  Cr  doped  10  (In2-xCrx03-5).  Part  1  of  this  section  will  include 
the  Sn  doped  10  system  and  will  provide  details  of  the  fabrication  methods,  structural  and 
optical  characterization.  Then  it  will  describe  an  extraction  methodology  based  on 
Nilsson’s  approach  [22]  on  these  films.  Part  2  will  present  our  most  recent  results  on 
magnetic  semiconductor  (Cr  doped  10)  based  homogeneous  negative  index  materials. 
This  section  will  provide  details  on  the  fabrication,  structural,  optical  and  magnetic 
characterization  of  In2_xCrx03_5  films,  and  provide  an  alternative  method  for  extracting  the 
negative  index  of  refraction  of  these  materials  based  on  reflection  and  transmission 
coefficients  measured  at  multiple  angles  of  incidence. 

1.  Sn  doped  IO  material  system 

This  section  describes  the  structural  and  materials  physics  of  the  doped  semiconductor 
(Ini_xSnx)203.  The  report  is  divided  up  into  three  sections  1)  a  description  of  the 
fabrication  of  the  (Ini_xSnx/2Znx/2)203  films  and  their  structural,  electronic  and 
compositional  characterization,  2)  a  report  of  the  results  of  electromagnetic  properties  of 
the  materials  in  the  ultra-violet/visible  spectral  range  and  the  mid-infra-red  and  3)  The 
computation  of  the  real  and  complex  parts  of  the  refractive  index.  The  real  part  of  the 
refractive  index  can  be  pushed  into  the  negative  by  pump  and  support  laser  fields, 
whereby  electrical  dipole  transitions  are  magnetically  coupled  to  produce  a  overlapping 
negative  regions  of  dielectric  permittivity  and  magnetic  permeability.  This  produces  the 
negative  index  effect  which  we  look  for  with  a  probe  laser  which  optically  samples  the 
region  of  the  film,  wherein  the  negative  refractive  index  is  produced. 

The  films  are,  at  best,  homogeneous  polycrystalline  smooth  slabs  of  material  with  carrier 
densities  and  conductivities  similar  to  those  seen  observed  in  n-doped  semiconductors. 
We  focus  on  a  set  of  a  series  of  In203  (10)  and  In2_xSnx03  (ITO)  films,  which  serve  as 
the  basis  set  of  samples  for  the  growth  and  development  of  the  In2-x(Sn,Zn)x03.  The  films 
are  characterized  by  scanning  electron  beam  (SEM)  microscopy  for  homogeneity  and 
smoothness,  surface  profiling  for  film  thickness,  energy  dispersive  X-ray  spectroscopy 
(EDX)  for  sample  composition,  Hall  measurements  for  electron  concentration,  optical 
transmittance  by  the  spectrometry  in  the  UV/ Visible  spectral  region  and  the  mid-infrared 
(mid-IR)  spectral  region  using  a  Fourier  transform  infrared  interferometer  (FTIR). 

The  optical  transmittance  data  can  be  fitted  to  a  Drude  plasmon  model  from  which  the 
high  frequency  permittivity  can  be  extracted,  along  with  the  plasma  wavelength,  the  bulk 
plasmon  damping  and,  in  principle,  the  effective  mass.  The  transmittance  and  reflectance 
data  are  then  used  for  computational  extraction  of  the  optical  constants  -  n  and  k.  All  thin 
film  samples  were  grown  by  sputtering  on  <100>  p-Si  substrates. 


1.1  Crystal  Structure  and  Film  Fabrication 

The  quality  of  the  sample  material  and  the  reproducibility  of  the  observed  optical 
properties  are  critically  important  to  the  success  of  this  project.  Over  the  last  three 
months,  a  significant  portion  of  time  has  been  devoted  to  the  development  of  the  material 
sources.  The  doped  (Sn,  and  Zn  in  this  instance)  ImCE  films  in  a  DC  and  RF  sputtering 
system.  Several  films  of  both  10  and  ITO  (Sn  doped  M2O3)  have  been  grown  with 
varying  Sn  concentration  and  Ar  gas  pressure.  The  Ar  gas  pressure  is  used  to  sputter 
from  the  target  and  we  have  found  that  a  pressure  of  24  mT  is  conducive  to  the  smooth 
and  homogeneous  films  whose  thickness  lies  in  the  range  of  several  microns.  This  is  ideal 
for  the  optical  transmittance  and  EDX  measurement  considerations. 
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Fig.  2  (a)The  cubic  bixbyite  structure  of  (II, VI)  randomly  doped  ln203 

The  thicknesses  are  measured  by  surface  profiling  and,  more  accurately,  by  measuring 
the  separation  in  adjacent  Fabry-Perot  (FB)  interference  fringes  in  the  transmittance  data. 
These  FB  oscillations  can  be  observed  in  the  mid-IR  range  for  films  of  several  microns  in 
thickness.  SEM  imagining  reveals  the  surface  morphology  and  smoothness  of  the  films. 

The  Sn  and  Zn  dopant  atoms  will  be  independently  controlled  in  terms  of  sputtering  rates 
and  hence  the  deposition  of  these  atoms  into  the  insulating  M2O3  matrix  is  done  in  order 
to  controllably  balance  the  density  of  these  two  atomic  species. 


Fig.  2  Sputtering  chamber  photograph  showing  the  sputtering  guns.  The  calibration  curve  enables  the 
control  of  the  composition  of  the  film  with  Sn/Zn  gun  power.  The  plots 
Show  the  growth  rate  calibration  for  deposition  rate  Vs.  power. 


1.2  Electrical  and  Compositional  Characteristics 

The  optical  scattering  and  losses  are  critically  dependent  upon  the  carrier  density.  This 
we  measure  using  a  Hall  setup,  from  which  the  Hall  constant  and  mobility  are  derived. 
The  typical  Hall  carrier  concentrations  that  are  measured  in  the  intrinsic  I112O3  are  in  the 
range  of  1x10  -5x10  cm'  .  This  concentration  is  pushed  upwards  by  the  addition  of  Sn 
atoms  into  the  material,  as  shown  in  Table  1. 


In  order  to  gain  an  understanding  of  the  composition  of  the  10  and  ITO  films,  we  subject 
our  films  to  EDX  analysis  with  a  field  emission  SEM.  This  beam  penetrates  a  distance  of 
about  1  micron  into  the  sample.  The  accumulated  statistics  on  collected  keV  dispersive 
energies,  reveal  the  Sn  (Z=49)  peak,  which  have  had  to  have  been  separated  from  the  In 
(Z=50)  peak.  From  these  atomic  percentages,  we  calculated  the  mole  ratios,  x,  and  the 
resulting  carrier  concentrations  as  measured  by  the  Hall  measurements.  We  find  that,  in 
general,  the  carrier  density  is  determined  completely  by  the  Sn  mole  ratio  (x)  and  the 
oxygen  mole  ratio  3-6  in  the  formula  In2.xSnx03+(5 . 


Table  1.  Summary  of  the  EDX  material  composition  analysis  of  the  DC/RF  co-sputtered  ITO  films.  The 


carrier  concentration  is  measured  by  the  Hall  effect. 


ln2-xSnx03 

sample 

In  at.% 

Sn  at  % 

X 

Ax 

density  (cnV-S) 

Oat  % 

thickness  him) 

sp10-016 

2643 

023 

0 01725431 

0  014381161 

9  61E+19 

45  86 

15 

sp10-017 

26  11 

1  36 

0  09901711 

0  014090541 

3  69E+18 

43  95 

1.17 

sp10-018 

61  05 

2.53 

0  07958478 

0  006066175 

5.84E+20 

20  26 

1.06 

sp10-019 

26  46 

0  97 

0.07072548 

0  014041874 

2.44E+21 

45.2 

0  619 

spl  0-020 

60  89 

3  04 

0  09510402 

0  006049774 

6 19E+20 

19  78 

1 

sp10-021 

42  37 

777 

0  30993219 

0  008530097 

9  80E+20 

16.89 

0  316 

spl  0-022 

67  79 

2  67 

0  07578768 

0  005470576 

3  41E+20 

22  95 

1  82 

spl  0-023 

60  54 

0  9 

0  02929688 

0  006243564 

1  20E+21 

21  49 

1  03 

1.3.  Optical  characterization  in  the  UV/visible 

In  order  to  understand  the  optical  characteristics  of  the  10  and  ITO  films  in  the  visible 
region  of  the  electromagnetic  spectrum,  we  performed  transmittance  measurements  on 
the  samples  using  a  broadband  halogen  and  deuterium  lamp  and  an  Avantes  3000 
visible/UV  spectrometer.  We  were  able  to  resolve  clear  Fabry-Perot  oscillations  in  the  IO 
films,  indicating  their  smooth  and  optically  homogeneous  quality.  The  UV/Visible 
transmittance  measurements  also  revealed  strong  and  pronounced  inter-band  absorption 
at  around  350nm  (see  Figure  3  (a)),  which  is  corresponds  to  a  band  gap  of  3.1-3.55  eV 
(Figure  3  (b))  and  the  absorption  occurs  between  the  between  O  (2p)  and  In  (5s)  states  in 
the  material. 
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Fig.  3  (a)  Transmittance  curves  of  the  indium  oxide  films  showing  Fabry-Perot  oscillations  from  which  the 
thicknesses  were  extracted.  The  dips  at  350nm  indicate  interband  absorption,  (b)  Optical  transmittance 
curves  for  the  ITO  films  -  showing  the  suppression  of  the  interband  absorption  due  to  scattering  losses,  (c) 
the  band  gap  estimate  linear  fit  for  a  typical  sputtered  IO  film  . 


Figure  3(b)  shows  the  plots  of  the  transmittance  for  the  ITO  films.  These  plots 
demonstrate  that  these  films  are  more  lossy,  due  to  the  inclusion  of  the  Sn  substitutions. 
The  interband  transitions  are  screened  out,  yet  the  Fabry-Perot  oscillations  are  still 
present,  indicating  that  the  films  are  smooth  on  the  length  scales  of  the  probing 
wavelength.  These  films  of  ITO  are  now  ready  for  the  doping  of  the  Zn  atoms.  The  Zn 


atoms,  unlike  the  Sn  atoms,  which  donate  one  electron  to  the  conduction  band,  are 
acceptors  and  give  two  holes  to  the  conduction  band.  The  Zinc  atoms  are  to  be  used  to 
reduce  the  total  carrier  concentration  while  maintaining  a  Sn  doping  mole  ratio  of  about 
x=0.05. 

1.4  Optical  characterization  in  the  mid-IR 

We  performed  transmission  measurements  on  the  10  and  ITO  films  using  broadband  THz 
radiation  with  a  Brueker  Optics  Vertex  70  FTIR  incorporating  a  room  temperature  Csl 
DLATGS  detector,  whose  sensitivity  falls  in  the  wavelength  range  2.5  /j m  to  30  /j m.  The 
substrate  material  p-Si  offers  a  50  %  transmittance  in  this  range.  This  substrate  is  used  to 
characterize  our  background  signal.  \  The  optical  constants  and  the  parameters  of  the 
electric  permittivity  function,  are  found  from  the  extraction  procedure  as  described  by  P. 
O.  Nilsson  using  the  formalisms  developed  by  Mayer  and  Heavens.  The  extraction  of  the 
Drude  parameters  are  done  by  fitting  the  transmittance  curves  to  a  Drude  model.  This 
latter  procedure  is  absolutely  necessary  to  establish  whether  we  have  negative 
permittivity  in  the  mid-IR  region. 

1.5  Drude  plasmon  fitting 

In  order  to  extract  the  theoretically  intractable  parameter  eM,  together  with  the  bulk 
plasma  wavelength  and  the  plasmonic  losses,  we  fit  the  mid-IR  transmittance  data  to  a 
model  of  the  permittivity  function.  This  is  done  by  writing  the  Drude  model  for  the 
dielectric  permittivity,  as  follows: 


Where  6cx3  is  the  high  frequency  dielectric  permittivity,  Ap  is  the  plasma  wavelength  and 
r  parameterize  the  losses.  The  transmittance  is  then  calculated  in  terms  of  the  Fresnel 
formulae  and  a  fit  with  the  mid-IR  transmittance  data  is  done.  A  typical  fit  is  depicted  in 
Figure  4,  while  the  table  of  fitted  parameters  are  given  in  Table  2. 


Fig.  4  A  typical  fit  of  the  mid-lR  transmittance  curve  for  an  10  sample  in  order  to  extract  the  3  Drude 
parameters  of  the  Drude  model.  The  plasmon  losses  are  underestimated  here  And  the  long  wavelength  fit  is 
divergent  from  the  data,  as  are  the  data  in  the  Fabry-Perot  interference  region. 


Table  2.  Drude  model  fits  to  the  transmittance  data  of  a  series  of  indium  oxide  films. 


Sample# 

einf 

XP(um) 

r 

Fit  converqence  (chi  sq.) 

09-012 

0 031413 

9  7854 

0.53037 

0  174118 

sp09-10 

1  294 

12  6477 

0  014763 

0  03 

sp09-13 

0  051615 

13  1977 

0  46964 

0  0342691 

sp09-14 

0  55895 

5  0104 

0  06071 

15.4424 

sp10-26 

0.051433 

14  7177 

0.5387 

0.0664292 

1.6  Numerical  extraction  of  the  optical  constants 

In  order  to  accurately  extract  the  optical  constants  of  the  films  (the  real  and  imaginary 
parts  of  the  refractive  index),  we  model  the  optical  transmittance  and  reflectance  using 
the  methods  described  in  Refs  22-24.  The  idealized  film/substrate/air  system  as 
considered  in  the  paper  can  be  represented  in  the  following  way: 


Air 


Air 


T  =  nsTf  ,  . 

/  I -0  -  r 


*  y 


The  equations  which  relate  the  film  properties  and  the  measured  transmittance  and 
reflectance  are  strongly  non-linear  and  implicitly  specified  as  follows: 

F(T,n,k)  =  0  and  G(R,n,k)  =  0 

The  right  hand  side  of  these  equations  is  taken  to  be  squared  and  this  is  the  functional  to 
be  minimized  for  the  extraction  of  the  optical  constants  of  the  film  in  question. 

The  Nilsson  extraction  procedure  is  very  powerful  for  the  ITZO  optically  pumped 
negative  index  experiment  because  all  of  the  optical  constants  of  a  film  can  be  extracted 
from  a  single  shot,  low  noise  transmittance  experiment  over  a  finite  spectral  range.  The 
phase  can  also  be  estimated,  eliminating  the  need  for  the  use  of  the  Kramers-Kronig 
relations,  altogether.  This  method  of  the  determination  of  the  optical  constants  for  the 
film  is  crucial  as  the  probe  laser  will  be  collinear  and  coincident  with  the  pump  and 
support  beams,  and  will  be  critical  for  accurately  determining  n  and  k  under  the  requisite 
Photon  pumping  conditions. 


2.  Cr  doped  IO  material  system 


In  this  section  we  will  discuss  in  detail  the  homogeneous  NIMs  based  on  magnetic 
semiconductor  In2_xCrx03.  We  have  shown  that  a  strongly  pronounced  negative  refractive 
index  effect  at  ~  27.8  pm.  This  effect  was  theoretically  predicted  earlier,  and  it  is  based 
on  coexistence  of  the  spin  wave  (magnon)  mode  with  the  plasmonic  mode,  with 
simultaneous  negative  permittivity  and  permeability  responses.  The  thin  films  of  Im.xCr 
XC>3  are  fabricated,  with  low  stoichiometric  oxygen  deficiency,  which  is  required  for 
ferromagnetic  behavior.  The  magnetic  measurements  clearly  reveal  the  ferromagnetism 
with  high  saturation  magnetization.  The  complex  frequency-dependent  refractive  index  is 
extracted  from  combined  transmittance  and  reflectance  data. 

Theoretical  studies  of  Cr  doped  IO,  [25]  predict  that  such  principally  homogeneous 
materials  (monocrystals  or  polycrystals)  belong  to  class  of  magnetic  semiconductors  with 
ferromagnetic  behavior,  and  detailed  calculations  of  the  expected  negative  refractive 
index  were  done  for  the  transparent  conducting  Cr-doped  indium  oxide  [26].  The 
calculations  show  that  this  ferromagnetic  system,  in  its  polycrystalline  form,  should 
possess  the  fully  isotropic,  low-loss  negative  refraction  index  n  <  0  at  ~  10.48  THz.  The 
negative  refraction  index  was  owing  to  overlapping  the  negative  permittivity  function 
e{cd)  <  0,  due  to  the  plasmons  (resonance  frequency  cop),  with  the  negative  permeability 
function  ju(co)  <  0  due  to  the  spin  waves.  In  support  of  the  idea  of  homogeneous  materials 
with  desirable  effect,  it  was  demonstrated  experimentally  [27],  that  ferromagnetic  metal 
mangenite,  La2/3Cai/3Mn03,  reveals  negative  refractive  index  at  -  150  GHz,  based  on 
analogous  plasmon-magnon  excitation  mechanism. 

2,1  Fabrication  of  Cr  doped  IO  material  system 

Fabrication  of ferromagnetic  In2-xCrx03-j  films  is  a  considerable  challenge  itself  since  the 
ferromagnetic  behavior  of  In2-xCrx03-5  critically  depends  on  the  stoichiometric  oxygen 
deficiency  concentration,  or  parameter,  8  which  should  be  within  the  difficult-to-reach 
narrow  range:  1.6x10"*  <  S<102  [28].  Since  every  oxygen  deficiency  provides  2  electrons 
into  the  conduction  band,  the  8  controls  the  electron  density,  ne,  in  conduction  band.  As 
the  ferromagnetic  ordering  in  this  material  is  electron-mediated,  and  heavily  depends  on 
ne_,  the  required  8  within  the  aforementioned  narrow  range  is  crucial  for  indirect  spin 
coupling,  which  causes  the  ferromagnetism  [28].  Also,  additional  details  of  the 
fabrication  process  dramatically  affect  the  magnetic  behavior  of  the  Cr-doped  IO,  and  for 
this  reason,  the  ferromagnetism  becomes  elusive  or  difficult  to  reach.  Moreover,  several 
groups,  for  this  very  reason,  even  deny  the  very  existence  of  the  ferromagnetic  properties 
[29-30]  or  suggesting  the  presence  of  secondary  phase  as  the  origin  of  ferromagnetism 
[31-32], 

To  cope  with  this  difficulty,  we  have  used  an  effective  post  annealing  process  for  precise 
control  of  structural  parameter  8.  The  material  growth  technique  used  for  depositing  the 
thin  films  was  based  on  the  DC/RF  magnetron  sputtering  technique.  The  sputtered 
samples  are  deposited  in  a  chamber  whose  base  pressure  was  maintained  at  around 
2xl0~7  Torr  and  the  films  were  deposited  on  p-Si  (100),  whose  resistivity  was  on  the 


order  of  0.1  Q  cm.  In  terms  of  substrate  temperature,  the  grown  samples  are  synthesized 
at  room  temperature. 


The  sputtered  samples  are  grown  using  two  targets  -  a  d.c.  magnetron  sputtering  gun  for 
the  Cr  and  an  r.f.  magnetron  gun  for  the  In2  O3  Cr /  In?  O3  targets  (K.J.  Lesker,  both  at 
purity  of  99.99%).  The  sputtering  is  performed  in  an  Ar  atmosphere  of  around  24  mT, 
while  the  oxygen  concentration  can  be  enhanced  by  the  inclusion  of  a  small  partial 
pressure  of  oxygen  at  around  0.3  x  10  3  mT.  The  oxygenation  comes  from  a  rf  driven  O2 
plasma  source  and  is  controlled  by  a  needle-valve  and  mass  flow  controller  5.8  seem.  In 
this  way,  we  have  grown  Cr  doped  indium  oxide  with  the  general  formula  of  In2  x  Crx03.g 
(ICO),  where  x  is  the  Cr  mole  fraction  and  8  represents  the  deviation  from  stoichiometric 
oxygen  content.  These  samples  were  post  annealed  at  250°  C  in  air.  Even  a  tiny  change  in 
the  wt%  of  the  oxygen  results  in  a  significant  change  in  the  oxygen  deficiency,  8. 

The  stoichiometric  composition  required  for  ferromagnetism  is  within  the  following 
range:  38.80±0.506  wt  %  In,  60.00±0.108  wt%  O  with  1.20±0.398  wt%  Cr.  In  order  to 
obtain  this  composition,  the  fabricated  samples  are  post  annealed  at  250  °C.  It  is  known 
that  annealing  in  vacuum  results  in  the  formation  of  oxygen  vacancies  [33-34],  On  the 
contrary,  annealing  of  CIO  in  air  causes  oxygen  to  re-  diffuse  into  the  film  and  eradicates 
the  oxygen  vacancies.  EDS  results  are  presented  in  Table  1  below.  As  follows  from  this 
table,  by  changing  the  time  of  annealing  (60  sec),  we  were  able  to  reach  the 
aforementioned  range  [28]  for  8. 

TABLE  I.  The  effect  of  post  annealing  time  on  the  composition  of  the  sputtered  Cr-doped  InO 


Annealing 

%0 

%Cr 

%In 

X 

6 

As-deposited 

75.91 

1.23 

22.85 

0.102 

-3.3048 

250°C  in  45  sec. 

59.40 

0.75 

39.85 

0.037 

0.0739 

250°C  in  60  sec. 

59.93 

0.73 

39.34 

0.036 

0.0087 

250°C  in  50  sec. 

64.14 

1.01 

34.86 

0.056 

-0.5762 

250°C  in  55  sec. 

62.89 

1.18 

35.93 

0.064 

-0.3894 

2.2  Characterization  of  Cr  doped  IO  material  system 

The  ordinary  Hall  Effect  with  a  Van  der  Pauw  contact  geometry  was  used  to  determine 
the  conduction  carrier  concentration,  ne.  The  structural  study  of  the  fabricated  samples 
was  done  by  energy  dispersive  spectroscopy,  EDS,  and  the  film  thickness  was  measured 
by  the  AlphaStep  profilometer  instrument. 

The  magnetic  measurements  of  the  successful  samples  are  performed  by  SQUID  at  10  K 
and  at  room  temperature,  in  hope  to  see  expected  ferromagnetic  behavior.  The  magnetic 
field  is  applied  in  the  sample  plane.  Fig.  1  shows  the  results  of  the  magnetization  versus 
applied  magnetic  field  (M-H)  curve  measured  at  10  K  for  the  0.35  pm  thick  of  Cr:InO 


annealed  thin  films  with  carrier  density  ne=  0.529  xlO19.  As  follows  from  this  result,  our 
successful  samples  are  indeed  ferromagnetic,  with  high  saturation  magnetization  up  to 
0.6  ps/Cr-atoms.  The  measurements  of  the  ferromagnetic  thin  samples  with  thickness  0.1 
pm  provided  even  higher  saturation  magnetization  up  to  0.8  pB/Cr-atoms.  Since  this 
saturation  magnetization  is  only  ~  25%  larger  than  magnetization  0.6  ps/Cr-atoms  for 
thick  samples,  we  can  conclude  that  magnetization  is  mostly  bulk  effect  and  surface 
enhances  ferromagnetism  only  to  a  small  extent. 
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Fig.  1.  Magnetization  measurements  indicating  that  the  doping  of  Cr  atoms  to  the  ln203  film  and  after  post 
annealing  yields  a  saturated  magnetic  moment  of  (0.6-0. 8)  pB  per  Cr  atom  at  10  K. 

Next,  the  transmittance  of  the  films  were  measured  at  normal  incidence  using  a  Bruker 
Optics  Vertex  80  interferometer  incorporating  a  broadband  THz  lamp  and  a  room 
temperature  KBr/DTGS-D301  photodetector.  This  system  enabled  us  to  perform 
transmittance  measurements  of  the  films  over  the  spectral  range  up  to  40  pm.  The 
reflectance  has  been  measured  within  broad  range  of  angle  of  incidence,  0 

(5,10,15,20,40,50,70,80  degrees)  by  the  AutoSeagull  reflection  unit. 

Fig.  2(a)  demonstrates,  for  a  ferromagnetic  sample,  the  examples  of  measured  frequency- 
dependent  FT1R  both  the  reflectance  coefficient  R  for  9  =  5° ,  (a)  0  =  50°  (b),  and  the 
transmittance  coefficient.  One  can  see,  from  Fig.2,  the  strongly  pronounced  maximum  on 
R  curves  at  4max-27.8  pm,  and  the  minimum  on  the  T  curve  at  the  same  frequency.  We 
assume  that  these  extremes  correspond  to  the  limiting  spin  wave  frequency,  oo  -  10.8 
THz,  on  the  boundary  of  the  Brillouin  zone.  Indeed,  this  frequency  is  very  close  to  the 
theoretical  predictions  [26].  Hence,  the  maximum  at  -  27.8  pm  is  a  viable  candidate  for 
the  desirable  effect,  and  we  extracted  the  refractive  index  within  narrow  range  which 
includes  this  wavelength  (see  below). 

It  would  be  instructive  to  compare  reflectance  spectra  of  ferromagnetic  and  non-magnetic 
films,  with  the  same  Cr  doping  parameter  x,  but  with  different  oxygen  deficiency,  8.  Fig. 
3  demonstrates  the  typical  FTIR  R  spectrum  for  the  non-magnetic  film,  with  the  same 
doping  x  -0.036  as  the  ferromagnetic  film  (8  -  0.0087)  utilized  in  Fig.  2,  but  with  the 
oxygen  deficiency,  8  ~  0.06,  outside  of  the  required  range  (1.6x10^  <8  <KT2)  which  is 
required  for  the  ferromagnetic  indirect  spin-spin  coupling  [28].  One  can  see,  in  contrast 


Field(Oe) 


to  ferromagnetic  film,  that  the  non-magnetic  film  possesses  no  maximum  in  region  27-28 
pm.  Such  behavior  is  expected,  since  the  spin  waves  (magnons)  which  are  responsible  for 
the  maximum  are  not  presented  in  this  specific  filmO.  Moreover,  the  non-magnetic  R 
spectrum,  at  27.8  pm,  shows  a  minimum  instead  of  a  maximum,  and  R  is  much  smaller 
than  R  for  the  ferromagnetic  film.  One  can  conclude  that  the  strongly  pronounced  peaks 
on  R  curve,  with  simultaneous  extremely  high  reflectance  R  >  0.8,  are  the  signatures  of 
the  magnon-plasmon  resonance  in  magnetic  semiconductors. 


R 
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Fig.  2.  Reflectance  coefficient,  R,  and  transmittance  coefficient,  T,  of  the  ferromagnetic  thin  film  of  In2-X 
Crx  O3.5,  x-0.036,  8  =  0.0087  ,  for  different  angles  of  incidence:  a  )  R,  9=5  ,  b)  R,  #=50  ,  c)  T,  9=0  ,. 
Red  arrow  points  to  magnon-plasmon  overlapping  resonance  at  ~  27.8  pm,  which  corresponds  to  a  maxima 

of  R  and  a  minima  of  T. 


R 


Fig.  3.  Reflectance  coefficient,  R,  of  the  nonmagnetic  thin  films  of  In2-xCrx  O3.5,  x~0.036,  S~  0.06,  for 

angle  of  incidence  #=50  . 


2.4  Extraction  Methodology 


Even  though  FTIR  measurements  do  not  provide  the  phase  of  reflected  and  transmitted 
waves,  the  optical  constants,  or  real  and  imaginary  parts,  n\(A)  and  112(A),  of  a  complex 
refractive  index,  11(A)  =  ii\(A)  +  i  112(A),  can  still  be  reliably  extracted  from  the  combined 
reflection  and  transmission  amplitudes,  R(A,6)  and  T(A,0),  measured  at  sequence  of 
incident  angles  0X ,O2,...0l  ,  if  this  sequence  covers  a  wide  range  of  angles  [36].  The 
method  of  the  extraction  of  n(A )  which  was  used  for  this  project,  is  based  on  the 
following  formulations.  Obviously,  each  theoretically  calculated  curve  with  fixed  0  for 
reflection  and  transmission  RTH(0,nx(A),n2(A)),  TTH(0,n1(A),n2(A)) ,  should  depend  on 

n\(A)  and  112(A).  Hence,  within  a  narrow  region,  close  to  the  specific  wavelength  of  an 
interest,  A~  Aext,  it  is  possible  to  fit  calculated  amplitudes 
Rm(6, nl(A),n2(A)),  Tm(6  ,nx(A),n2(A))  to  experimental  data,  R(A,6)  and  T(A,0),  in  order 

to  extract  the  couples  of  optical  constants  {nl(6i,A),n2(0i,A)}  which  depend  on  the 
incident  angle  6=  0  x.  Afterwards,  one  can  analyze  which  couple  [nx,  n2 }  is  the  same  for 
all  R-T  amplitudes,  or  for  the  whole  set  6  =  6n  and  this  solution  should  correspond  to 
true  refractive  index,  n^(Aext),  n2(Aexl )  . 

In  the  calculation  of  the  reflection-transmission  amplitudes,  Rth(&,  «i (A),  112(A))  and 
Tth(6,  n\(A),  112(A)),  we  have  followed  the  classical  Bom’s  method  [36]  of  characteristic 
reflection-transmission  matrix  of  a  stratified  medium.  This  medium  is  the  film  which 
covers  the  silicon  carbide  substrate  with  thickness  -  250  pm  and  permittivity  es  ~  3.38  + 
/0.0034  (both  extrapolated  from  data  of  Ref.  37  and  directly  measured  at  wavelength  -28 
pm). 

The  appropriate  extracted  refractive  index  band  in  vicinity  of  the  magnon-plasmon 
resonance,  at  Amax  -  27.8  pm,  is  shown  in  Fig.  4.  One  can  see,  from  Fig.  4,  that  refractive 
index  becomes  negative  within  narrow  band  with  Re(n)  ~  -2.0,  \m(n) «  2.0  ,  which 
corresponds  to  figure  of  merit,  FOM  -  1 .  Hence,  due  to  negative  refractive  index  effect  at 
~  27.8  pm,  we  can  conclude  that  this  wavelength  consistently  corresponds  to  assumed 
plasmon-magnon  resonance. 

As  follows  from  Fig.4,  the  positive  refractive  index,  Re(«)  ~  1.0  and  \m(n)  ~  10.0, 
outside  of  the  negative  refractive  index  band.  The  magnitude  of  this  refractive  index  can 
be  explained  from  the  Drude  plasmon  pennittivity  s  [22], 

s(co)  =  £j\  -  0)0 >- — -)  (1), 

(co  +iyco) 

where  the  Drude  frequency,  cor>~  310.5  THz,  corresponds  to  concentration  of  electrons 
(N~2.7x  1 0"  (1/cm  ))  in  the  conduction  band  for  experimental  oxygen  deficiency  8  ~ 
0.0087;  ~0.8,  and  yis  the  losses  in  plasmonic  subsystem.  Due  to  the  lack  of  magnetic 

response  (//  -  1)  outside  the  narrow  band,  at  the  boundary  of  the  Brillouin  zone  [25,  26], 


(which  coincides  approximately  with  the  negative  refractive  band),  the  refractive  index, 
outside  the  band,  can  be  approximated  as  n  =  \[s  .  Since  the  Drude  frequency  is  much 
larger  than  the  frequency  where  the  negative  refractive  index  band  is  located  (ojd 
»  cam),  the  calculated  real  part  of  the  Drude  permittivity  is  negative:  Re(£)  — 100.0. 
The  appropriate  real  and  imaginary  parts  of  the  index  of  refraction  are  as  follows:  R e(n)  ~ 
1.0  and  Im(n)  ~  10.0.  Hence,  as  one  can  see  from  Fig.  4,  real  and  imaginary  parts  of 
refractive  index,  as  predicted  by  Drude  theory,  are  on  the  same  order  of  magnitude  with 
experimental  refractive  index  in  the  vicinity  of  the  negative  refractive  index  band.  The 
behavior  of  n(A)  is  fully  consistent  both  within  and  outside  the  negative  refractive  index 
band  with  the  theoretical  predictions  [26,  38]. 


Fig.  4  Extracted  refractive  index  «=Re(«)+/Im(«)  in  the  vicinity  of  the  plasmon-magnon  resonance  of  the 
ferromagnetic  film  In2-xCrx  03_s,  x-0.036,  S~  0.0087. 


In  conclusion,  the  negative  refractive  index  band  parameters  of  ferromagnetic  Cr-doped 
indium  oxide  thin  films  (i.e.  wavelength  and  bandwidth)  are  fully  consistent  with 
predicted  ones  in  Ref.  27,  and  the  measured  electric  and  magnetic  properties  of  Cr-doped 
IO  are  close  to  those  reported  in  literature  [29].  Specifically,  the  experimentally  verified 
negative  refractive  index  narrow  band,  with  width,  AAJa lmax~0.005,  is  located  at  the 
limiting  frequency  of  the  magnon  spectra,  2max~27.8  pm,  with  Rc(«)~  -2.0,  Im(/?)  -2.0. 
The  theoretical  prediction  of  Ref.  10  provides  the  negative  refractive  index  band  at  Amax~ 
30.0  pm  with  the  estimated  width,  H2,//lmax~0.00 1-0.1,  and  the  refractive  index,  Re(«)~  - 
2.5,  Im(ft)  -1.0.  These  predictions  are  in  close  proximity  to  the  extracted  experimental 
parameters  reported  above. 


Conclusions 


As  a  result  of  this  project,  we  have  investigated  two  different  systems  for  achieving 
homogeneous  negative  index  materials.  Firstly,  doped  semiconductors  exposed  to  special 
external  laser  beams  to  achieve  coherent  coupling  of  an  electric  dipole  transition  with  a 
magnetic  dipole  transition  leading  to  permeability  and  pennittivity  responses  which  result 
in  negative  index  in  a  particular  wavelength  regime  is  considered.  Then,  magnetic 
semiconductors  (specifically,  Cr  doped  10)  are  considered  where  the  effect  is  due  to  the 
coexistence  of  the  spin  wave  mode  with  the  plasmonic  mode.  Both  of  these  modes  are 
activated  by  the  e.m.  field  of  the  light  with  simultaneous  permittivity  and  penneability 
responses  within  some  frequency  band,  which  ensures  the  negative  refractive  index 
within  the  frequency  band  close  to  the  boundary  of  the  Brillouin  zone  of  the  magnon 
spectra.  Based  on  these  studies  presented  herein,  we  believe  that  natural  homogeneous 
magnetic  semiconductors  with  well-pronounced  negative  refractive  index  band  can  be 
promising  in  future  applications.  The  advantages  of  these  natural  materials  compared 
with  inhomogeneous  composite  metamaterials  are  their  optical  isotropy,  and  the  fact  that 
the  optical  constants,  s,ju,n  are  true  physical  variables  defined  on  the  atomic  level, 
rather  than  some  “effective”  parameters. 
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